Vibrational Coherences in Nano-Elastic Tunneling 
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Charging a nano-scale oscillator by single electron tunneling leads to an effective double-well 
potential due to image charges. We combine exact numerical diagonalizations with generalized 
Master equations and show that the resulting quantum tunneling of the mechanical degree of freedom 
can be visualized in the electronic current noise spectrum. 
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Nanoelectromechanical systems (NEMS) display a va- 
riety of interesting dynamical effects due to the coupling 
between mechanical and electronic degrees of freedom. 
Some recent examples are oscillators with linear coupling 
to single [1] or multiple Q electronic levels, or systems 
with non-linear coupling such as quantum shuttles [3]. 

A common feature of NEMS is the possibility to study 
quantum- mechanics in the true meaning of the word, i.e. 
quantum coherent behavior of objects (molecules or me- 
chanical resonators |4( ) that are big on the scale of smaller 
units (electrons, atoms) [Hj]. In this Letter, we show that 
coherent quantum tunneling occurs in a small vibrating 
resonator that is driven by stochastic image charges and 
a corresponding switching between harmonic and non- 
harmonic (Coulomb-like) potentials generated by single 
electron transport. For certain parameters, the resonator 
tunnels back and forth between two spatially separated 
positions which can be visualized by measuring the noise 
spectrum of the electrons. 

Model .- Our model resonator oscillates in x-direction 
with frequency uoq and mean square displacement Zq m 
a harmonic oscillator potential that is modified by tun- 
neling of a single electron through a nearby quantum 
dot or molecular level. Coulomb attraction of the res- 
onator towards an image charge induced at a distance 
2|a| in a nearby surface occurs by charging with single 
electrons (cf. FigO]), e.g., in a molecular junction, in 
quantum dots embedded into a vibrating beam [6] , or in 
larger quantum dots where charging with an additional 
electron attracts the electron 'droplet' towards a nearby 
metallic gate. We thus assume two electronic 'dot' states 
|0) and |1) (with zero or one transport electron), and an 
additional Coulomb potential on the resonator in the oc- 
cupied state, leading to the effective Hamiltonian of the 
dot-resonator 
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where M = H/IqUUq is the effective oscillator mass of the 
resonator, u the unit vector from the dot to the image 
charge, and d) the creation operator of an electron on 
the dot, where for simplicity we disregard the electron's 
spin degree of freedom here and in the following. 




FIG. 1: Quantum dot in a nanoscale resonator close to a 
surface: single electrons (de) charging a nanoscale resonator 
at rates Tl(Xr) l ea d to a switching between a harmonic and 
a non-linear potential due to the Coulomb attraction of the 
charged resonator towards a surface image charge. 



We start our calculation by an exact diagonalization 
of iJdot • I n the three-dimensional problem, the resonator 
motion in x-direction decouples from the motion perpen- 
dicular to the x-axis for both charged and uncharged 
states, thus leaving us with an effectively one-dimensional 
problem: the uncharged vibrational states |n, 0) simply 
are harmonic oscillator eigenstates with respect to the 
equilibrium position a. Non-harmonic corrections to the 
uncharged resonator potential in Eq. (pQ) can easily be 
incorporated but do not lead to drastic modifications of 
our results. 

On the other hand, the charged states |n, 1) are ob- 
tained by expansion into the eigenbasis of I = hydro- 
gen atom eigenfunctions (weighted Laguerre polynomi- 
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sequent numerical diagonalization. Here, x is measured 
from the surface in units of a. 

The model has two dimensionless parameters, 
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where Q determines the distance from the wall relative 
to the harmonic oscillator (zero-point fluctuation) length 
scale Zq? an d e is given by Q times the typical Coulomb 
(charge - image charge) energy, 
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relative to the harmonic oscillator phonon energy HcjUq. 
Both parameters determine the shape of V(x) (cf. Fig. 
Efa)) and the resulting vibrational eigenstates |n, z/), 
where v — \(y — 0) denotes (un)charged resonator states 
and n is the vibrational quantum number. In the follow- 
ing, we measure all energies (rates) in units of (h)coo/Q. 

An important feature of the corresponding spectrum of 
the energy eigenvalues e vn for v = 1 is the appearance of 
an avoided level crossing (Fig. [2fa) inset near ft = 13.4). 
This can be traced back to the effective double- well struc- 
ture of V(x) for coupling strengths e > 9. In this regime, 
the charged oscillator can be approximated by a two-level 
system at low energies, and one thus expects quantum co- 
herent behavior due to tunneling between two states ap- 
proximately localized around x = a and x = (i.e. near 
the surface). In the following, we demonstrate that this 
vibrational tunneling can be made visible by electronic 
transport through the resonator, using Wigner function 
representations of the stationary density matrix and the 



frequency-dependent electronic current noise spectrum 
S(lo). 

Method .- We use the usual tunneling Hamiltonian 
Ht = J2ka Vk,ac{d + H.c. with left and right (a = 
L,R) electronic reservoirs at chemical potentials pl 
and hr described by the free electron Hamiltonian 
#res = J2ka £ k c {a c ka' We describe quantum transport 
by formulating a generalized Master equation for the 
reduced dot-resonator density operator p t that is suit- 
able for strong Coulomb blockade (two electronic states) 
and weak reservoir coupling. This approach describes 
only first order tunneling processes, i.e. neglecting co- 
tunneling effects. Importantly and in contrast to the 
usual rate equation method [l|, [7[ , we fully retain all vi- 
brational coherences in the resulting Master equation, 
which in the usual Born-Markov approximation is ob- 
tained by taking vibrational matrix elements of the Li- 
ouvillian [8[, leading to 
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and a corresponding equation for the uncharged states. 
The indices j and f run over the uncharged oscillator 
states and i over the charged states. Note that no co- 
herence between charged and uncharged states occurs as 
must be. Here, we defined the rates j a (s) = T a f a (e) for 
tunneling into the dot, the rates 7 a (e) = r a (l — f a (^)) 
for tunneling out of the dot where the coupling strength 
to the leads is defined as T a = Vj? a 27r5(e — Also, 
we use a short hand notation for the matrix elements 
p nu m^ = (n,v\ pt \m, p) and f a (e) for the Fermi function 
at p a + £• 

The Franck- Condon factors naturally occur by 
sandwiching p t with the vibrational states. They are de- 
fined as the overlap integral of the vibrational wavefunc- 
tions (x \n,v) = (j) n ,v{x) before and after the transition: 

nj = / dx (j)i^(x)(j)j^(x). (5) 
Jo 

The Tij play an important role for strong vibrational cou- 
pling e. Here and in contrast to the usual small polaron 
models (linear oscillator coupling) p], 0, 0] , they have to 
be calculated numerically from the exact eigenfunctions 
of i^dot- Some Franck-Condon factors are shown in Fig. 
Efb). As a function of £7, there are two different regimes 
where either m or T21 is dominating. The sudden rise of 
Tn is due to the level anti-crossing after which the shape 
of the first eigenstate resembles the harmonic oscillator 



ground state and vice versa for T21 and the second eigen- 
state. Also the sign of the has to be taken into account 
(cf. Eq. [4j) contrary to the usual rate equation approach 
where the Franck-Condon factors occur as squared val- 
ues only. Our numerical evaluation shows that the vi- 
brational off-diagonal elements of Eq.(|4j) are not small in 
the regimes we are interested in and therefore must not 
be neglected. 

Results .- Even without damping of the resonator (see 
below), our system has quite a large parameter space 
(£1, e, pr, r^, Tr) and consequently exhibits a rich be- 
haviour. We therefore first visualize the stationary state 
Pstat = p(t — > 00) of the dot-resonator system by plot- 
ting its Wigner transform, W(x,p) = ^- dye zpy (x + 
^y\p\x — \y), which displays the expected behaviour as 
mentioned above: the system clearly 'lives' in two sepa- 
rate regions in phase space, one being directly near the 
surface (x = 0) and well localized in x direction, whereas 
the other region strongly depends on the parameters but 
is generally located around the equilibrium point a (Fig. 
[3]) . The concentric rings around this point are the contri- 
bution of the empty oscillator, which are just the usual 
harmonic oscillator Wigner functions. These contribu- 
tions 'interfere' with the contributions of the charged os- 
cillator which, as a result, leads to the two- level charac- 
teristics visible in the total Wigner function. 

The stationary current, Fig. [2](c)&(d) shows the usual 
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FIG. 2: (color online), (a) Potential of an oscillator with 
image charge effect depending on the oscillator frequency Q 
and the eigenvalue spectrum of the correspondong Hamil- 
tonian (for e = 12) (inset). Level anti-crossing occurs at 
Q = 13.4. (b) Franck-Condon overlaps of the eigenfunc- 
tions of the charged and uncharged oscillator as a function 
of Q. There are two different regimes before and after the 
anti-crossing point, (c) dl/dV as a function of gate and bias 
voltage for e = 12 and H = 11 without the effect of the 
Franck-Condon factors (nj = 1) showing the usual conduc- 
tance diamonds, (d) As (c) but with Franck-Condon factors: 
Current is surpressed for low bias and cannot be restored by 
changing the gate (Franck-Condon blockade). Negative dif- 
ferential conductance occurs as violet lines in the transport 
spectrum (c.f color scale) 



signatures of a Franck-Condon blockade, i.e. a suppres- 
sion of the low bias current that cannot be lifted by 
application of a gate voltage and which is due to the 
Franck-Condon factors that describe the overlap of 
vibrational states before and after tunneling of electrons. 
Furthermore, we also observe negative differential con- 
ductance [10] (Fig. [2(d)), and huge Fano factors [l| (Fig. 
Ha)). 

The most significant information about the resonator, 
however, can be directly extracted from the electronic 
noise spectrum S(u). Following the method developed 



U|, we calculate the frequency dependent noise from 



the convenient expression 
Srr(u) = (1 R ) 
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where Xr is the part of the Liouvillian that generates a 
quantum jump (an electron tunneling to the, e.g., right 
reservoir), 1Z{uo) is the pseudo-inverse of the Liouvillian 
C and (•) = Tr[m 

Pstat]- Note that for simplicity, we 
only calculate the noise contribution in the right reser- 
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FIG. 3: (color online). Wigner distribution of the dot- 
resonator for coupling strength e = 12 and different parame- 
ters Q, cf. Eq. J2]). The contribution of the harmonic oscil- 
lator Wigner distribution (uncharged resonator states \n, 0)) 
is clearly visible as a series of concentric rings around the 
equilibrium point at x = 1. For small Q, the system is well 
localized close to the surface, whereas for intermediate Q it 
clearly lives in two separate regions of phase space correspond- 
ing to the two effectively localized two- level states. For larger 
Q, there is again localization around x — 1. 



voir which (together with the contribution from the left 
reservoir) is sufficient to obtain the noise for strongly 
asymmetric capacitance coefficients. 

Figs. H(b)&(c) demonstrate that the noise is indeed 
extremely sensitive [12| to changes in the energy level 
spectrum of the resonator. The level splitting close to 
the avoided level crossing becomes visible as the appear- 
ance of a feature in S(co) at frequencies uo close to the 
corresponding Bohr frequencies (level splitting in the vi- 
brational spectrum, Fig. [4(b) inset). The electronic noise 
thus acts as a direct probe of vibrational coherence. We 
point out that the appearance of non-trivial spectral in- 
formation in S{uj) is due to retaining the non-diagonal 
elements in the Master equation Eq. (jlj) . In order to em- 
phasize this point, we have compared our results with a 
corresponding calculation where these elements were set 
to zero (Fig. H(b)). 

In order to elucidate the role of damping, we again 
sandwiched the usual Lindblad form £damp = \ (2apa^ — 
pa) a — a)ap) with the exact vibrational eigenfunctions. 
Here, 7 is the damping rate of the uncharged oscillator, 
described by a and a) . The most dramatic effect this 
damping has on the system is a considerable reduction 
of the zero frequency Fano factor (Fig. HI (a)). The huge 
zero frequency Fano factor can be understood as a dy- 
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FIG. 4: (color online), (a) Zero Frequency Fano factor for 
zero damping (lower frame) and several damping strength 
7 £ [10 -3 , 8 x 10 -3 ] (upper frame), (b) Frequency dependent 
Fano factor F(u) = S{uo)/(X) for systems with eigenfrequen- 
cies Q as indicated in the inset, VBias = VGate = 30. Coupling 
parameters are Yl = 1/2tt and Fr = l/4n. For comparison 
the case where only diagonal elements of the Liovillian are 
considered (black line), (c) As (b). The sensitivity of the 
noise to the level anti-crossing is seen as a growing and flat- 
tening peak at the frequency corresponding to the level spac- 
ing for eigenfrequencies scanning through the anti-crossing 
point. The oscillator eigenfrequency of each line is indicated 
by the right frame as well as its position with respect to the 
avoided- crossing. The scale of the axis refers to the oscilla- 
tor frequency at the level anti-crossing point = 13.4 (red 
line) the other lines have an offset of =bn0.05. (d) F(ui) for 
£1 = 13.4 and several damping strength 7. For high damp- 
ing the peak at the level spacing disappears and the noise 
spectrum resembles the single resonant level result. 



namical channel blockade as described in [l|: The peaks 
observed in Fig. HJa) all occur at Q values where at least 
one Franck-Condon factor is zero (cf. Fig. [2](b)). This 
means one transport channel is blocked and the current 
is withheld for some time until it can flow through some 
other channel. The peak at ft = 9.6 is unusually high 
compared with the other because here several Franck- 
Condon factors are zero or of the order of 1(T 2 . The 
damping now has the effect of allowing relaxations be- 
tween the levels and thus even for very small damping 
the channel blockade is lifted. 

The observed resonance peak at the level spacing ener- 
gies is also suppressed for a system subjected to damping 
(Fig. Ufa)), and the quantum coherence between the vi- 
brational states is gradually lost. For high damping con- 
stants 7 the whole noise spectrum approaches the limit 
of a static quantum dot with two levels [13[ , as one would 
expect. 

Finally, we comment on possible experimental realiza- 
tions of our predictions. The two model parameters Q 
and e, Eq. ([2]), fix the ratio M/cj = h 3 (87re s r /e 2 ) 2 e 2 /n 



between effective oscillator mass M and frequency ujq. 
The values Vt = 13.4 and e = 12 used above yield a dis- 
tance a ~ 10 _9 m and would in fact correspond to an 
'electron droplet' oscillator of a few electron masses, i.e. 
quantum dot electrons collectively oscillating at frequen- 
cies uoq ~ 10 14 Hz (using e r ~ 10). In calculations for 
larger e, however, we also confirmed that one can easily 
reach M of the order of a few proton masses or larger, 
with the two-level approximation still being valid. These 
masses are reasonable for vibrating molecular junctions 
or even larger macroscopic oscillators. Since Q is pro- 
portional to the distance a it could be tuned by using a 
metallic STM tip as the part of the environment were the 
image charge occurs. Then one would be able to scan Q 
and thus find the two level parameter regime, where the 
described noise features can be observed. 
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